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Using first principles calculations, we investigate the electronic and magnetic properties of the 
trigonally distorted (i?3c) perovskite-derived NiCrOa. Within the local spin density approximation 
(LSDA), our calculations show that this system is an exactly compensated half- metal (CHM). The 
local spin moments of Cr 2.04, and antialigned Ni -1.41 and three oxygens -0.63 (in the units of 
/is), indicate high spin S — Cr^^ and 5* = | (NiOs)'^" units. Considering reasonable values of the 
on-site Coulomb repulsion U on both Ni and Cr ions with LDA-I-U approach, this system becomes 
an insulator (as reported by Chamberland and Cloud) having a narrow gap in the spin-up channel, 
whereas the other channel has a large gap of ~3 eV. Although inclusion of U seemingly leads to 
the transition Ni^""" — >-high spin 5* = | Ni"^"*", consistent with the experimentally observed effective 
moment, the zero net moment remains unchanged due to either reduction of oxygen local moments 
or enhancement of Cr local moment. Compression of volume by 10% leads to CHM even when 
correlation effects are included. These results suggest the possibility of a CHM state in NiCrOa and 
provide another route to search for CHM, which is a property sought by many. 

PACS numbers: 71.20.Be, 75.47.Lx, 71.27.+a 



Introduction. A half-metallic antiferromagnet, 
more properly called a compensated half-metal 
(CHM), is half- metallic, having one conducting spin 
channel and one insulating spin channel, but van- 
ishing (compensating) net magnetic moment, hence 
no macroscopic magnetic field. Since proposed by 
van Leuken and de Groot[l| fifteen years ago, the 
CHM has been considered as an encouraging candi- 
date for spintronics - fully polarized carriers with- 




FIG. 1: (Color online) Crystal structure of the trigo- 
nally distorted NiCrOa, having the corner-sharing CrOe 
octahedra. 



out a macroscopic magnetic field to deal with - and 
a possible single spin superconductor. [Ij However, a 
true CHM has not yet been established. 

Since at least two different kinds of magnetic ions 
are necessary for vanishing moment, the double per- 
ovskite systems were investigated in several theoret- 
ical investigations. [3-0] However, preparin g a n or- 
dered sample has been frequently failed. [Tol.TlTI| The 
Heusler structures, which have produced many half- 
metals, also have been considered as another promis- 
ing candidate. [T2l - [l3 | Nakao has suggested CHM in 
tetrahedrally coordinated chalcopyrites[ll| and in 
their monolayer super lattices. [l6j 

In this study, we will extend the search for CHMs 
to the trigonally distorted perovskite-related sys- 
tems having two different transition metals. Several 
systems of this class were synthesized about forty 
years ago. [T7l - [l9j Here we will investigate NiCrOa 
with a corundum-related structure, synthesized by 
Chamberland and Cloud, [l^l Their procedure in- 
volved preparation at high (60-65 kbar) pressure at 
1200 °C, followed by quenching to room tempera- 
ture, so intermixing on the Ni and Cr sites cannot be 
ruled out. The susceptibility measurements showed 
a Curie- Weiss behavior in the high T regime and 
discontinuities at 250 K and 120 K, proposed to be 
canted antiferromagnetism. The value of the effec- 
tive moment was Pe//=5.8 ^b, close to the value 
for Cr^+ and Ni'^^, both carrying high spin 5=1, 
implying a possibility of CHM. This system was pro- 
posed to be semiconducting with an activation en- 
ergy of 0.11 eV, though no detailed information was 



provided. More recently, this compound has been 
obtained as a minority during the oxidization pro- 
cess by a few groups, [13, who measured only the 
space group and lattice constants consistent with 
those of Chambcrland and Cloud. [l7| 

Structure and methods. In this trigonally dis- 
torted system, Chambcrland and Cloud reported 
only the space group i?3c (No. 167) with the hexag- 
onal lattice constants a = 4.9252 A and c = 13.504 
A,[l3 leading to a distortion factor of 0.89 

(unity for the ideal perovskite). We have optimized 
the structure parameters in the cubic and the trig- 
onally distorted structures. In both structures, the 
one in which the Ni ion sits on the A site and the Cr 
ion on the B site is much favored over reversed site 
occupation. Consistent with the experiments, our 
optimizations show that the cubic phase with opti- 
mized lattice constant a = 3.602 A has much higher 
energy, by a few eV, than the distorted phase. Hence 
we will consider only the distorted structure. 

In the i?3c structure, using the experimental lat- 
tice parameters, [l7j our optimizations indicate that 
Ni, Cr, and O ions lie on the 6a (0,0,i), 66 (0,0,0), 
and 18e {xfl,j) sites, respectively, with the oxygen 
internal parameter x of 0.3748. As displayed in Fig. 
m the Cr-0 bond length of 1.915 A is about 3% 
larger than the Ni-0 bond length of 1.85 A. The 
CrOe octahcdra consist of a nearly square base and 
equilateral triangular faces, with only 1.8% differ- 
ence between the 0-0 bond lengths. Thus the CrOe 
unit remains close to an ideal octahedron. The O- 
0-0 bond angles are 60°, 2x60.6°, and 58.8°. The 
0-Cr-O bond angles are either 91° or 89°. The Ni 
ions lie at the center of mass of the equilateral trian- 
gle framed by in-plane oxygens. This trigonal pla- 
nar symmetry D^h of the NiOa leads to the crystal 
field splittings: doublet eig {dyz,dzx), doublet e2g 
{dxy,dx2-y2), and singlet aig {d^^). 

These calculations were carried out within the 
local spin density approximation (LSDA) followed 
by the LDA-I-U approach, as implemented in the 
all-electron full-potential code FPLO-9.[22 In the 
LDA+U approach, two popular double-counting 
schemes, the so-called around mean field (AMF) and 
the fully localized limit (FLL),[2l ^ were used, 
showing similar results for this compound. In our 
calculations, a rhombohedral cell containing two for- 
mula units was used. The convergence was checked 
with a regular mesh containing up to 781 irreducible 
k points. 

LSDA electronic structure. The Cr local moment 
(2.04 fis) is precisely canceled by the moments of 
antialigned Ni (-1.41 /is) and three oxygens (-0.63 
/is). The oxygen moment is unusually large for a 
3d system, though not unprecedented. (In the PBE 
generalized gradient approximation, [2^ the CHM 





FIG. 2: (Color online) Top: LSDA band structure near 
the Fermi energy Ef, where the bands arise only from 
d orbitals of Ni and Cr ions. The dashed horizontal line 
indicates Ef, which is taken as the zero of energy. Bot- 
tom: Brillouin zone and high symmetry points of the 
rhombohedral structure, which has similarity to the fee 
Brillouin zone with a threeefold axis oriented along the 
z axis. 



configuration also results but with somewhat larger 
moments: Cr 2.26 /xb, Ni -1.54 ^b, three oxygens 
totaling to -0.72 fis ■) This CHM state with zero 
net moment is favored by 300 meV/f.u. over the 
nonmagnetic state. Thus NiCrOs is a compensated 
half-metal within LSDA, [2^ regardless of the choice 
of exchange-correlation functional. 

The LSDA band structure near Ep is displayed in 
Fig. [2l illustrating the 0.8 eV gap in the spin-down 
channel leading to the half-metallic state. In the 
spin-up (metallic) channel, the Ni e2g bands lying in 
the range of 0.6-1 eV are separated from the mix- 
ture of the other Ni d and Cr t2g bands spread out 
over the regime of -1.8 to 0.4 eV. In the spin-down 
bands, the Ni 3d (majority) states are fully occupied 
with small (1.4 eV) bandwidth, and Cr 3d (minor- 
ity) states arc empty. The corresponding densities of 
states in the full energy range and with enlargements 
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FIG. 3: (Color online) Total and atom-projected LSDA 
densities of states, which shows the half-metallic aspect, 
expressed per eV- formula unit. In all panels, the shaded 
regions indicate three oxygens-projected DOS. The total 
density of states N{Ef) at Ef is 3.80 states per eV per 
a formula unit, decomposed of 50% Ni, 30% Cr, and 20% 
three oxygens characters, (a) Cr 3d DOS in solid line, 
from -2 eV to 2 eV; (b) Ni 3d, on same scale as in (a); (c) 
total DOS in the entire O 2p - metal 3d region. The un- 
occupied Cr Cg states with the exchange splitting of 1.5 
eV lie on 1.5-4 eV and 2.5-5 eV in the spin-up and spin- 
down channels, respectively. The exchange splitting of 
Cr t2g states is about 2 eV. The t2g-eg crystal field split- 
ting of Cr in the spin-up channel is 3.5 eV, 1 eV larger 
than in the spin-down channel. The vertical dashed line 
indicates Ef- 



emphasizing Ni and Cr separately, are displayed in 
Fig. 131 The oxygen bands in the range of -7.5 to -2 
eV are separated from the d bands by 0.5 eV. How- 
ever, the mixture of O 2p character into both Ni and 
Cr bands indicates strong p — d hybridization. 

The spin-up channel is more interesting, being a 
mixture of Cr majority and Ni minority states. To 
balance the five occupied Ni (majority) spin-down 
bands (per f.u.) there must be five spin-up bands 
occupied (on average through the zone, due to the 
metallic nature). Three bands, around 0.5 eV and 
above, are separated from the other intermixed seven 
bands, so the filling is 5/7. These can be identified 
as Cr Cg and a substantial admixture of Ni (and of 
O 2p as noted earlier). The occupied spin- up bands 
are Cr t2g admixed with Ni dy^, dzx, and bands. 
Being strongly hybridized and metallic, formal va- 
lences become a murky concept. S = ^ Cr'^+ seems 
to be the clearest feature; presuming ions with 
their formal valence one would obtain •S' = | Ni'^+. 
We emphasize that the strong p — d hybridization 
and metallicity make these assignments of limited 
use. 

Inclusion of correlation effects. We have used the 
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FIG. 4: (Color online) Blowup band structures along 
the Ux-X-F line near Ef for various values of U in both 
schemes of LDA-I-U. Here, n6c4 denotes U=6 and 4 eV 
for Ni and Cr ions, respectively. In panel (b), the gap of 
5 meV separates nearly degenerate linear bands similar 
to what has been observed (at the F point) in CoSbs.fS^ 




E-Ep (eV) 

FIG. 5: (Color online) Total and atom-projected den- 
sities of states per formula unit in the AMF scheme of 
LSDA+U using [/= 4 and 6 eV for Cr and Ni ions, re- 
spectively. Near Ef the occupied part has mostly Cr d 
character, but the bottom of the unoccupied part has Ni 
d character even though the strong O characters also ex- 
ist. In particular, the O characters are nearly the same 
as the d character in the range of —0.1- 0.1 eV. This 
indicates a nearly compensated half-semimetal. 

LDA+U approach to model the observed insulat- 
ing state. The proper U values for this system are 
unclear, but values that have been widely used in 
oxides are 3-4 eV for Crjl^ and 6-8 eV for Ni,!! 
having been obtained from comparison with experi- 
mental data or calculated from the constrained LDA 
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approach. [30| The Hund's exchange integral J=l eV 
has been fixed for all values of U, after some varia- 
tion of J indicated results are insensitive to values 
around this one. 

Figure |4] displays the enlarged band structures 
along the Ux-X-T line (which is where the gap 
opens) near Ep (set to zero), for these ranges of 
U and using both commonly used LDA+U schemes. 
The energy gap in the spin-down channel increases 
to about 2.6 eV in AMF and 3.6 eV in FLL (not 
shown in the figure). In the spin- up channel, the 
bottom of the conduction bands and the top of the 
valence bands both consist of the mixture of the O 
Pz and Ni e2g bands, in particular near the X point, 
with little Cr d character. However, about 50% at 
the bottom of the conduction bands and the top of 
the valence bands are oxygen and Ni characters, re- 
spectively. The second highest valence band, which 
is nearly flat, is one of the Cr t2g bands. In AMF, 
at n6c4 (denoting Uni = & eV, Ucr = 4 eV) a small 
gap in the spin-up channel opens at the X point, 
whereas in FLL two bands almost touch ai Ep- For 
n8c3, in FLL the band structure is similar with that 
of nQcA, except for a slight increase in the energy 
gap. In AMF, however, Cr t2g bands at the F point 
are raised, leading to an small indirect gap of the F 
to X point almost equal to the direct gap at X. This 
indirect gap is observed only in the case of U eV 
for Cr ions in AMF. 

The FLL and AMF methods are known to en- 
courage magnetism differently. [3l| Here we find that 
applying U on the metal atoms leads to significantly 
different induced moments on the oxygen ions. In 
FLL, applying U enhances the magnitudes of all lo- 
cal moments, and the combined moment of three 
oxygens is still about 40% of the Ni moment. In 
AMF, on the other hand, both Ni and Cr moments 
are increased by applying [/, but the moment of 
three oxygens becomes only 20% of the Ni moment. 
These different moments can be ascribed to an en- 
ergy penalty for magnetism in AMF that is not 
present in FLL, as discussed previously. [sij On bal- 
ance, however, the moment compensation remains 
unchanged when varying U in the range studied 
here. 

Since differences in the overall DOS curves are 
small, we show only that for 7i6c4 in AMF in Fig. 



O In the spin-up channel, the Fermi energy lies be- 
tween the occupied Cr d bands and the unoccupied 
Ni e2g bands. Just above 1.5 eV, narrow unoccupied 
Ni aig bands, i.e. mixing slightly with O states, ex- 
ist. As a result, in the insulating phase, the Ni ion is 
apparently trivalent, even though the local moment 
of -1.90 /is is considerably reduced due to strong 
p—d hybridization. However, the charge difference of 
Ni ions, obtained by the Mullikan charge decompo- 
sition, between in LSDA and LDA-I-U is only ~0.2e, 
similar to that observed in charge disproportionated 
Nai/sCoOs. 
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Discussion and Summary. Given the very small 
gap in the up channel, modest pressure should close 
the gap, giving a CHM state, as we now show. At 
10% compression in volume, the zero moment is 
slightly degraded within LSDA, but inclusion of U 
in the same range as in the ambient pressure re- 
stores the state of the zero moment. In the spin- up 
channel, which is insulating at ambient pressure, Ni 
aig bands of W «0.6 eV cross Ep along the L-V-X 
line, resulting in being metallic and a precise com- 
pensated half-metal. (Below ^ 10% compression, 
the electronic structures are very similar with those 
at ambient pressure.) This suggests that applying 
pressure will lead to a CHM state even though this 
system is insulating at ambient pressure. 

In summary, we have investigated i?3c-structure 
NiCrOs; using both LSDA and LDA-I-U approaches. 
Within LSDA, this system is a precise compensated 
half-metal with the energy gap of 0.8 eV in the spin- 
down channel. Considering the on-site Coulomb re- 
pulsion U of both Ni and Cr ions, we obtained an 
insulating phase, but retaining zero net moment re- 
gardless of the magnitude of U . Interestingly, the 
energy gap is substantial (2.6 - 3.6 eV) in the spin- 
down channel, whereas it is tiny or vanishing in the 
spin-up channel. Application of pressure provides a 
likely route to a CHM state in this compound. These 
findings provide another structure class that may be 
favorable for synthesis of a compensated half metal. 
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